clear
rng(20)

simulated_data_baseline = [];
simulated_data_low_spillover = [];
outer_counter = 1;
repeats = 1;
num_comps = 10;
%parpool('local',num_comps)

%inputs = [0.9; 0];

%%% BASELINE PARAMETERS
%%% EXPORT_SLOPE = 2.35
params = [0.1207	0.0296	0.1584	1.538	7.108	0.938];
inputs = [export_slope; 0];
iteration_data = zeros(1,7);
counter=1;
while counter <= repeats
parfor computer_num = 1:num_comps
temp_data = steady_state_simulation_export_creation_destruction(params, inputs);
iteration_data = iteration_data+temp_data;
end
counter = counter+1;
end
iteration_data = iteration_data./(repeats*num_comps);
simulated_data_baseline = [simulated_data_baseline; export_slope iteration_data]

export_creation_destruction_rate_kappa_094 = iteration_data(1)

%%% LOW SPILLOVERS
%%% THETA AND ARRIVALS RE-ESTIMATED
params = [0.1653 0.0344 0.2467 1.538 6.9405 0.01];
export_slope = 1.82
inputs = [export_slope; 0];
iteration_data = zeros(1,7);
counter=1;
while counter <= repeats
parfor computer_num = 1:num_comps
temp_data = steady_state_simulation_export_creation_destruction(params, inputs);
iteration_data = iteration_data+temp_data;
end
counter = counter+1;
end
iteration_data = iteration_data./(repeats*num_comps);
simulated_data_baseline = [simulated_data_baseline; export_slope iteration_data]

export_creation_destruction_rate_kappa_001 = iteration_data(1)

